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I ■ ABSTRACT 

o 

■^ , Context. With the latest infrared surveys, the number of massive protostellar candidates has increased significantly. New studies have 

t^ . posed additional questions on important issues about the formation, evolution, and other phenomena related to them. Complementary 

to infrared data, radio observations are a good tool to study the nature of these objects, and to diagnose the formation stage. 

Aims. Here we study the far-infrared source IRAS 16353^636 with the aim of understanding its nature and origin. In particular, we 

search for young stellar objects (YSOs), possible outflow structure, and the presence of non-thermal emission. 
►^ . Methods. Using high-resolution, multi-wavelength radio continuum data obtained with the Australia Telescope Compact Array* , we 

QQ ' image IRAS 16353-4636 and its environment from 1.4 to 19.6 GHz, and derive the distribution of the spectral index at maximum 

/«^ ' angular resolution. We also present new JHK^ photometry and spectroscopy data obtained at ESO NTT**. '^^CO and archival Hi line 

p^ \ data, and infrared databases (MSX, GLIMPSE, MIPSGal) are also inspected. 

iy~, ■ Results. The radio continuum emission associated with IRAS 16353-4636 was found to be extended (~ 10 arcsec), with a bow-shaped 

• ' morphology above 4.8 GHz, and a strong peak persistent at all frequencies. The NIR photometry led us to identify ten near-IR sources 

r^ ' and classify them according to their color We used the Hi line data to derive the source distance, and analyzed the kinematical 

(«_^ I information from the CO and NIR lines detected. 

<«_^ . Conclusions. We have identified the source IRAS 16353-4636 as a new protostellar cluster. In this cluster we recognized three 

distinct sources: a low-mass YSO, a high-mass YSOs, and a mildly confined region of intense and non-thermal radio emission. We 

propose the latter corresponds to the terminal part of an outflow. 

Key words, infrared: stars - ISM: individual objects: IRAS 16353-4636 - radio continuum: stars - stars: pre-main-sequence 

1. Introduction in finding thousands of candidates. The goals of these surveys 

have been to confirm the evolutionary status of luminous, em- 
The availability, in recent years, of infrared survey products like bedded sources, to perform a statistical analysis of these objects 
the MSX (Price et al. 2001) and succesive catalogs have pro- ^^ g^j^^tj^ s^^jes^ ^n^j ^^ ^^^y. ^j efficient mechanisms for mas- 
vided an extremely rich database to look for massive young ^j^„ „f^ formation 
stellar objects. Programs such as the Red MSX Survey (RMS; 

Lumsden et al. 2002, Hoare et al. 2005) have been successful 

Observations of star forming regions and associated outflows 



* The Australia Telescope Compact Array is funded by the ^^ radio wavelengths can be used to further our understanding 

Commonwealth of Australia for operation as a National Facility by of massive star formation, a topic which is still not completely 

CSIRO. understood. However, such studies require high angular resolu- 

** Based on observations collected at the European Organisation tion and sensitivity, because massive young stars are often found 

for Astronomical Research in the Southern Hemisphere, Chile (ESO at kilo-parsec distances and are usuaUy associated with densely 

Programme 073.D-0339, PI S. Chaty). populated clusters of intermediate and low mass stars. 
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Up to now, only a few massive young stellar objects (massive 
YSOs or MYSOs) are related to collimated jets mapped at the ra- 
dio range, such as the Serpens sources (Rodriguez et al. 1989), 
HH 80-81 (Marti et al. 1993), and IRAS 16547-4247 (Garay et 
al. 2003). Garay et al. found this last object to be a triple quasi- 
linear radio source, that shows non-thermal indices at the lobes. 
The source could be associated with a molecular core. The au- 
thors have proposed that it is a MYSO ejecting a collimated wind 
that interacts with the surrounding interstellar medium, produc- 
ing shocks and consequent radio emission. In recent years, near- 
IR and (sub)millimeter line studies have proved very valuable 
in supplying more information on (M)YSOs; see, for example, 
Varricatt et al. (2010) and references therein. 

We have found an intense infrared source, 
IRAS 16353-4636 (l,b = 337.99°, 0.08°) which seems to 
share some characteristics with IRAS 16547-4247, the object 
observed by Garay et al. (2003). The source IRAS 16353-4636 
has been cataloged as a star-forming region by Avedisova 
(2002), based on its IRAS colors, and studied for the first time 
by Combi et al. (2004). Although the available multi- wave length 
data from radio to high-energy gamma rays suggested it could 
be a microquasar candidate, later X-ray observations with 
XMM-Newton provide a more precise position that was away 
from the IRAS source and its radio counterpart (Bodaghee et 
al. 2006). In fact, the X-ray source has been found to be an 
accreting X-ray pulsar 

The main goal of this study is to establish the nature of 
IRAS 16353-4636. Section 2 describes the already known 
sources in the target field, and details the observation and re- 
duction processes we have carried out. In Sect. 3 we present the 
results. A discussion is given in Sect. 4 and the conclusions are 
stated in Sect. 5. 



2. Observations and data reduction 

2.1. The field of IRAS 16353-4636 

The infrared source is centered at {a,S) [J2000] - 16*^39'" 
3.52^ -46°42'28.28". The uncertainty ellipse major and mi- 
nor axes and position angle are 14" x 4", 97° (Beichman 
et al. 1988). Mid-IR and Far-IR fluxes are given in Table 1. 
Sources from other catalogs are found in the neighborhood 
and we plotted them in Fig. 1. The MSX6C catalog (Egan et 
al. 2003) lists the source G337.9947h-00.0770, which coincides 
with IRAS 16353-4636 and is detected at the four MSX bands 
(Table 1). 

In the frame of the RMS Survey, Urquhart et al. (2007) have 
detected '^CO (J=l-2) emission lines from the MSX source, us- 
ing the Mopra Telescope (angular resolution: 20", velocity reso- 
lution: 0.2 km s ', and noise temperature: 0.1 K). Very recently, 
Mottram et al. (2010) have estimated the MIPSGal 70yum flux of 
the RMS source (Table 1). 

In addition, there are two Spitzer GLIMPSE point sources 
(SSTGLMC G337.9986-hO.0758 and G337.9907+0.0733, 
Benjamin et al. 2003, Churchwell et al. 2009) in the vicinity 
of our target, though neither of them is positionally coincident 
with IRAS 16353-4636. 

The X-ray source IGR J 16393 -4643 (Bodaghee et al. 2006, 
Chaty et al. 2008, Corbet et al. 2010), which is of X-ray pulsar 
origin, is also plotted in Fig. 1. The region also hosts a low- 
frequency radio source, J163903. 9-464215. 55, detected with 
the Giant Millimetre Radio Telescope (GMRT; Pandey et al. 
2006) at 610 MHz. There is no spatial correlation between 
IRAS 16353-4636 and the X-ray and GMRT-radio sources. 



Table 1. Relevant infrared sources in the field. 



Source 


a,6in000] 


A 


Flux 


Size/ 




(iims,dms) 


(fim) 


(Jy) 


ang.res 


IRAS 


16 39 03.5, 


12 


13.5+0.7 


14" X 4" 


16353-4636 


-46 42 28 


25 

60 

100 


80.5+4.0 

<806 

2210+330 


RA.=97° 


MSX 


16 39 03.4, 


8.28 


6.339+4.1% 


18" X 18" 


G337.9947 


-46 42 27 


12.13 


9.454+5% 




+0.0770 




14.65 
21.34 


8.788+6.1% 
35.59+6% 




MIPSGalt 


16 39 03.4, 

-46 42 27 


70 


698.24+5.24 


2" X 2" 



f: Flux value from Mottram et al. 2010. 

Table 2. ATCA observational parameters, and integrated flux 
density at each frequency. 



V 


Conf. 


^int 


Weight"' 


Synth, beam. 


rms 


^ int 


(GHz) 




(min) 




Position angle 


(mJy/b) 


(mJy) 


1.384 


6A 


247 


R=0 


7.'7x5'.'9, -3.3° 


1.00 


90+4 


2.368 


6A 


247 


R=0 


4'.'6x3'.'7, -4.3° 


0.50 


2+4 


4.800 


6A 


201 


R=0 


2'.'4xl'.'6, -16.6° 


0.14 


72+3 


8.640 


6A 


201 


R=0 


l'.'4xa'9, -19.2° 


0.09 


31+5 


17.344 


6C 


528 


R=2 


l'.'0xa'7, +1.4° 


0.10 


81+20 


19.648 


6C 


528 


R=2 


l'.'0xa'6, -0.2° 


0.09 


65+15 



Note: mJy/b: mJy beam ' . tt: R = Robust, see text. 
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Fig. 1. IRAS 16353-4636 (ellipse) and other sources in its 
field. Spitzer-GLIMPSE point sources: filled circles (Benjamin 
et al. 2003); MSX6C G337.9947+00.0770: square, GMRT-610 
MHz J163903.9-464215.55 (Pandey et al. 2006): hollow circle; 
IGR J 16393 -4643 (Bodaghee et al. 2006): larger cross. SmaUer 
cross: observations of '-^CO line (Urquhart et al. 2007). 



2.2. K[CA 

We have carried out radio continuum observations toward 
the source IRAS 16353-4636, with the Australia Telescope 
Compact Array (ATCA), at the frequencies 1.384, 2.368, 4.800, 
8.640, 17.344, and 19.648 GHz. 
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Fig. 2. ATCA radio continuum images of IRAS 16353-4636. A: Image at 1.4 GHz; contour levels of -3, 3, 5, 10, 18, 30, and 45 
times the rms noise of 1.0 mJy beam"'. B: Image at 2.4 GHz; contour levels of -3, 3, 5, 10, 18, 30, and 45 times the rms noise of 
0.5 mJy beam"'. C: Image at 4.8 GHz; contour levels of -3, 3, 6, 9, 12, 20, 30, 45 and 60 times the rms noise of 0.14 mJy beam"'. 
D: Image at 8.64 GHz; contour levels of -3, 3, 6, 9, 12, 20, 24 and 28 times the rms noise of 0.09 mJy beam"'. E: Image at 17.344 
GHz; contour levels of -3, 3, 5, 8, 11, 14, 17, 19, and 20 times therms noise of 0.10 mJy beam"'. F: Image at 19.648 GHz; contour 
levels of -3, 3, 5, 8, 11, 14, 17, 19, and 20 times the rms noise of 0.09 mJy beam"'. Synthesized beams are plotted at bottom left 
corners. North is up and east is to the left. 



Observations in tlie cm range. The data were taken in 2004 
January 11, in the 6A configuration, at full synthesis (16:00- 
4:00 UT). The field was observed interleaving simultaneous 
observations at two frequency pairs: 1.384/2.368 GHz and 
4.800/8.640 GHz. PKS 1934-638 served as the flux calibrator, 
and the nearby source 1646-50 was observed for phase calibra- 
tion. The total bandwidth at these four frequencies was 128 MHz 
over 32 channels. The total time on source at each frequency re- 



sulted in 247 min at 1.384 and 2.368 GHz, and 201 min at 4.800 
and 8.640 GHz. 



Observations in tine mm range. We observed simultaneously 
at 17.344 and 19.648 GHz in 2006 March 29, using ATCA in the 
6C configuration (13:00-0.30 UT). The sources PKS 1934-638, 
1253-055 and 1646-50 were used as flux, bandpass, and phase 
calibrators, respectively. The total bandwidth was 128 MHz over 
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32 channels. The total time on source was 528 min. The high 
phase stability allowed us to perform 12-min target source scans. 

All six data sets were reduced and analyzed with the miriad 
package. The calibrated visibilities were Fourier-transformed us- 
ing 'natural', 'uniform', and 'robust' weightings. The robust 
weighted images combine the 'natural' lower noise and the 'uni- 
form' higher angular resolution (Briggs 1995). The robustness 
parameter was set to R - Q for the 6A (cm) data, and to +2 for 
the 6C (mm) datcQ. The resulting synthesized beams and rms are 
given in Table |2] 

2.3. VLA data 

We searched for data in the VLA archives and found observa- 
tions centered relatively close to IRAS 16353-4636, taken in 
2001 February 1 at 20 cm in the BnA configuration. The data 
were taken with one IF centered at the rest frequency of the 
hyperfine transition of Hi (1420.406 MHz), with 127 channels 
over a total bandwidth of 1.56 GHz, giving a width of 12.2 kHz 
(or 2.58 km s ') per channel. After Hanning smoothing, this re- 
sulted in a velocity resolution of 5.16 km s ' per channel. The 
spectral observations were centered at vlsr = km s"' . The data 
were edited, calibrated, and imaged using the software package 
Astronomical Image Processing System (AIPS) of NRACQ. The 
resulting synthesized beam was 10" x 3'.'8. 



2.4. NTTNIR 

We conducted photometric (J, H and Ks filters) and spectro- 
scopic (0.9-2.5 yum) NIR observations of IRAS 16353-4636 on 
2004 July 10 with the spectro-imager SofI, installed on the ESO 
New Technology Telescope (NTT). The large field imaging of 
Sofl's detector was used, giving an image scale of 0'.'288 pixel"' 
and a field of view of 4:94x4^94. 

We repeated a set of photometric observations for 
each filter with 9 different 30" offset positions including 
IRAS 16353-4636, with an integration time of 90 seconds for 
each exposure, following the standard jitter procedure that al- 
lows a clean subtraction of the blank sky emission in NIR. 

The IRAF (Image Reduction and Analysis Facility package) 
suite was used to perform the data reduction, which included 
flat-fielding and NIR sky subtraction. For the three images, one 
in each filter, we obtained an astrometric solution by using more 
than 400 coincident 2MASS objects, with a final rms of O'.'O? in 
each coordinate. 

We carried out aperture photometry and transformed the in- 
strumental magnitudes into apparent magnitudes with the stan- 
dard relation: magapp - maginst - Zp - ext x AM, where magapp 
and maginst are respectively the apparent and instrumental mag- 
nitudes, Zp is the zero-point, ext the extinction, and AM the 
airmass. The observations were performed through an airmass 
close to 1 . 



' 'Robust' weights are a function of local u-v weight density: in re- 
gions where the weight is low [high], the effective weighting is natural 
[uniform]. The optimal value thus depends strongly on the configura- 
tion. Larger values of R can produce beams that have nearly the same 
point source sensitivity as the naturally weighted beam, but with en- 
hanced resolution. 

^ The National Radio Astronomy Observatory is a facility of the 
National Science Foundation operated under cooperative agreement by 
Associated Universities, Inc. 



The spectroscopic observation consisted of 12 spectra 
with the Blue and Red grisms of 1000 resolution, and a 
wavelength coverage of 9000 to 25000 A. The position of 
IRAS 16353-4636 in the slit was off^set 30" in half of the ex- 
posures to subtract the blank NIR sky. The total integration time 
was 240 s in each grism. We took Xe lamp exposures to perform 
the wavelength calibration. The 5'-long slit was centered at one 
of the most intense NTT sources that overlaps the radio emission 
{a,6[J2QQQ] = 16''39"3.18%-46°42'3L48"). The NTT source 
was labeled #5 (see Sect. 3.3 and Fig. 5). The position angle of 
the slit was -1-15° (positive from north to east), and its width, 1". 

The NIR spectra were reduced with IRAF by flat-fielding, 
correcting the geometrical distortion using the arc frame, shift- 
ing the individual images using the jitter off'sets, combining these 
images, and finally extracting the spectra. The analysis of SofI 
spectroscopic data, and more precisely the sky subtraction, was 
difficult due to a variable sky, mainly in the red part of the blue 
grism, causing some wave patterns. We observed a photometric 
standard star Hip084636 (G3 V) from the Catalog of Persson et 
al. (1998), to correct for telluric absorption. The spectra were 
finally shifted to the lsr rest frame. 

3. Results 

3.1. ATCA 

The continuum images built with the ATCA data are shown in 
Fig. 2. The angular resolutions range from ~ 8", at the low- 
est radio frequency, to less than 1" at the highest one. The 
radio source is detected at all frequencies. The position of the 
peak flux (PFP) at all frequencies agrees among position uncer- 
tainties. At the highest angular resolution (19.6 GHz) we mea- 
sure (a, (5)pfp[J2000] = 16'>39™3.53\ -46°42'29.3". The inte- 
grated flux density Sint is derived summing all the flux above 
the 3cr contour (Icr - image rms). We adopted an error in the 
integrated flux density equal to the absolute difference between 
the integrated flux above Icr and Smt- Table 2 lists the results. 
Throughout the text we have followed the convention Sv oc v'', 
i.e. a spectral index a > for thermal emission, and a < for 
non-thermal emission. 

At the lowest frequency -1.4 GHz, and the largest beam- no 
structure is appreciated. At 2.4 GHz extended emission is de- 
tected toward the N-W, though it is better defined at 4.8 GHz. 
The detections from 4.8 to 19.6 GHz show rather elongated 
sources, in the direction SW-NE. The images at 17.3 and 19.6 
GHz reveal a clumpy structure. 

The estimate of the spectral index can help to characterize 
the radiation regime and, eventually, can lead to a firm identifica- 
tion of the astrophysical source. We build a spectral index map at 
the highest angular resolution possible, i.e., using the data at 17.3 
and 19.6 Ghz. To attain the same angular resolution at the two 
frequencies we restore the images to the same Gaussian beam 
size (1" X 1") for both data sets. 

We consider only pixels with signal-to-noise > 5. In Fig. 3 
we present the resultant spectral index map. It displays fine-scale 
structure on size scales comparable to the NTT infrared sources, 
as shown by the K^ contours superimposed on to the spectral 
index map in Fig. 3. 

In general, the index is negative (non-thermal) towards the 
radio peak seen in Fig. 2, and towards the south of the infrared 
cluster associated with IRAS 16353-4636, and is positive (ther- 
mal) towards the north and west. 

Radio observations from massive YSOs have revealed both 
thermal (as in Marti et al. 1995) and non-thermal outflows or 



Benaglia et al.: The star forming region IRAS 16353-4636 



-46°43'22"-: 




16''39"04.00' 



03.75' 03.50' 03.35' 03.00' 

Right Ascension (J3000) 



Fig. 3. Spectral index map derived for the 
17.344 GHz/19.648 GHz ATCA data, in grayscale: black 
is -1, white is +1. The NTT /Tj-band emission is superimposed 
as blue contours of 25, 50, 100, 150, 200, 300, and 400 Jy. The 
grey contour corresponds to 3cr level continuum emission at 
17.344 GHz. NTT sources here named 5 and 8 are marked (see 
text). 




-100 100 

LSR Radial Velocity (km/s) 

Fig. 4. 2 1 -cm H i absorption spectrum for the target, plotted as a 
function of lsr radial velocity. The measured spectrum is divided 
by the continuum level and represents e^'^, where t is the optical 
depth. The horizontal dashed line is drawn at 5 iine/5 continuum - 1 
and the vertical dashed line is drawn at an lsr radial velocity of 
Okms'. 



jets (like in Rodriguez et al. 1989). Strong shocks can occur at 
the end points of the jets, giving rise to diffusive shock particle 
acceleration which, in turn, will produce non-thermal emission 
of synchrotron origin (see, for instance, Romero 2010, and ref- 
erences therein). 

The rough matching of the 17.3/19.6 GHz spectral index 
distribution with the NTT image in Fig. 3 shows that (i) the 
strongest of the southern sources (here called #5, see below) is 
superimposed on a radio emitting region with a « 0; (ii) at the 
position of the strongest of the eastern sources (#8, see below), 
a > 0; and (iii) a negative spectral index corresponds to the pfp 
radio maximum. 



3.2. VLA 

We produce a continuum image at 1 .42 GHz using the line-free 
channels in the +5 to -n 152 km s"' range. It basically matches the 
ATCA image obtained at 1.384 GHz, thus confirming the results 
of the ATCA observations albeit with pointing differences. These 
are probably due to poor phase referencing at low elevations for 
the VLA. 

From the VLA data set, we also obtain Hi spectra toward 
the peak in these radio data. The center of the box region over 
which the H i emission was averaged coincides with star #5 in 
our near-IR observations (see below). However, due to the larger 
angular resolution of the VLA data, the region encompassed all 
the near-IR sources relevant here. The H i spectrum of the target 
is presented in Fig. |4] Hanning-smoothed to a resolution of 5.2 
km s"'. At negative velocities there is broad absorption down 
to -100 km s"', followed by a detached feature spread over 4 
channels, centered at a velocity of approximately -120 km s ' . 

3.3. NTTNIR 

We show the final J, H and Ks band images in Fig. 5; we also 
compare our 17.3 GHz radio map with the Ks band image in 



this figure. Besides diffuse NIR emission, there are several point- 
like objects. We have identified the brightest ones, and numbered 
them from 1 to 9 at increasing right ascension (labeled in Fig. 5, 
bottom left panel). 

The daophot package in the IRAF suite has been used to ex- 
tract all the sources in the crowded field through PSF fitting. 
The aperture has been adapted from 4 to 10 pixels, up to the 
Airy's ray if possible, to each source, only to take into account 
the stellar flux, subtracting the nebula flux. At the end of the ex- 
traction process, we apply an aperture correction so as to finally 
use the same aperture photometry for all of the sources. An ad- 
jacent annulus of 1 pixel outer radius is used to estimate the sky 
background. The results are quoted in Table [3] The magnitudes 
of the whole extended emission (point-like sources and diffuse 
emission) are 7 = 13.1+0.5,//= 11.9 + 0.2 and /T^ = 10.7 + 0.2. 

The NTT photometry is used to build a color-color (C-C) 
plot. (J-H) versus (H-Ks) colors are plotted in Fig. 6. As a ref- 
erence, we have included the loci of dwaif , giant and supergiant 
stars, following the intrinsic colors listed by Tokunaga (2000). 
Straight dotted lines show the reddening vectors for these stars, 
and define a band where stars with normal colors, only affected 
by reddening, are expected to be found in the diagram. All the 
objects lying to the right of this band show an infrared excess in 
their H - K^ colors and are presumably young stars still un- 
dergoing some form of accretion. 

From a global point of view, all the objects in Fig. 6 seem 
to share a common amount of reddening, equivalent to 8 to 10 
magnitudes of extinction in the visual band. The colors of ob- 
jects #1, 2, 3, and 4 are characteristic of reddened normal stars 
like foreground objects. The rest of the sources shows an in- 
frared excess at K^-hand; they could represent an embedded clus- 
ter of YSOs. We have examined in particular the two brightest 
Kfi sources with colors of protostellar objects (stars #5 and #8). 
We use NTT spectroscopy, and SED fitting to study #5 and #8, 
respectively. Notably, source 8 lies close to the IRAS and MSX 
peak in this region, and is the most extreme outlier in the C-C di- 
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Fig. 5. Upper left panel: /-band image of the region toward IRAS 16353-4636, obtained with the NTT on 2004 July 10. Levels: 
1, 3, 5, 10, 20, ..., 100, 120, 150, 190, 250, 320, and 400 Jy. The uncertainty in the position of the IRAS source is indicated by an 
ellipse. Upper right panel: The same as the upper left panel but in the H band. Levels: 5, 10, 20, ..., 400 Jy. We show the position 
and size of the slit used for spectroscopic purposes. Lower left panel: The same as the upper left panel but in the ^s band. Levels: 
20, 30, ..., 400 Jy. The extended emission appears more clearly than in the J- and //-band images. Several point-like sources are 
visible, and labeled with numbers as in Table 3, at increasing right ascension. Lower right panel: Continuum emission at 17.344 
GHz in greyscale and grey contours, superimposed on to black contours corresponding to strongest near-IR /Tj sources. North is up 
and east is to the left. 



agram. Source 8 also coincides with a region of positive spectral 
index in our radio data (Fig. 3). This is consistent with thermal 
radio emission which could be associated with a thermal radio 
jet or optically thick HII region (Marti et al. 1993). Source 8 
is therefore potentially the best candidate for an intermediate or 
even high-mass YSO in the cluster. Source 5, on the other hand, 
coincides with non-thermal radio emission, as one might expect 
if this source was also a YSO associated with an outflow (see 
for example Bosch-Ramon et al. 2010 for theoretical considera- 
tions). It is unlikely, however, that source 5 is a UCHII region, 
since UCHII regions are usually associated with thermal emis- 
sion (Wood & Churchwell 1989). 

Figure 5, bottom right panel, shows the correlation of near- 
IR emission [K^ band) with radio emission at 17.3 GHz. The 
radio peak has no near-IR counterpart. The comparison tells 
us that the spectral index/NTT sources coincidences evident in 
Fig. 3 must be taken with caution, since the radio emission as- 
sociated with NIR sources 5 and 8 appears to be swamped by 
extended emission associated with the radio peak itself. Sources 
5 and 8 are not discretely resolved in the radio data. So its not 
entirely clear that the negative and positive spectral indices ob- 
served around sources 5 and 8 are really linked to these NIR 
point sources. The bulk of the radio emission may be unrelated 
to the NIR sources. 



3.4. Individual sources 

Source #5. Our NTT spectrum of source 5 reveals the presence 
of H I emission lines and absorption in the CO vibrational bands 
('^CO and '^CO, Av = 2 and 3), typical of pre-main sequence 
stars (see Nisini et al. 2005 and references therein). Figure 7 
shows the spectrum with individual lines identified in the J, //, 
and Ks bandfl 

We find the following spectral lines: Pa/?, Pa(5, '^CO lines and 
bands, '^CO bands, Bry, and possibly Br 10, Br 1 1. All of them 
are characteristic of star-forming regions. Many lines peak at ve- 
locities close to or more negative than the H i absorption feature 
at -120 km s"'. Table 4 lists the suggested line identifications, 
measured wavelengths, rest wavelengths, and measured lsr cen- 
tral velocities. A standard Galactic rotation circular model, such 
as Brand & Blitz (1993), assigns a minimum velocity of about 
-125 km s ' to the gas at the tangent point. Table 4 indicates 
tliat tlie velocity shifts are extremely widespread. 

Errors in the determination of radial velocities are strongly 
dependent on the signal of the individual lines measured and 
the wavelength resolution of the spectrum. For relatively strong 
emission lines, such as Bry, the uncertainty in the estimation of 
the centroid of the Gaussian profile is expected to be about 1/10 
of a pixel, which is equivalent to about 14 km s '. The uncer- 



We checked the lines at |http : //www . j ach . hawaii . edu/UKIRT/astrono 
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Table 3. NIR coordinates, the band at which they are measured (7-band: 1.25;um, //-band: 1.65/^m, /Ts-band: 2.17 jum), and magni- 
tudes of the 9 objects, ordered by increasing right ascension and marked in Fig. 5, lower panel. 



Star* 


Q'[J2000.0] 
(h,m,s) 


5[J2000.0] 

r,',") 


Band 


J 


H 


K, 


1 


16 39 02.84 


-46 42 26.99 


H 


15.12 + 0.04 


13.26 ±0.04 


12.41 ±0.07 


2 


16 39 02.87 


-46 42 28.09 


H 


15.19 ±0.05 


13.35 ±0.05 


12.56 ±0.09 


3 


16 39 03.02 


-46 42 24.49 


H 


17.11 ±0.11 


14.98 ± 0.06 


13.61 ±0.08 


4 


16 39 03.01 


-46 42 26.58 


K, 


16.26 ±0.07 


14.18 ±0.06 


13.37 ±0.14 


5 


16 39 03.18 


-46 42 31.47 


H 


15.28 ±0.05 


13.78 ±0.05 


12.63 ± 0.09 


6 


16 39 03.25 


-46 42 32.92 


J 


15.30 ±0.05 


13.93 ±0.05 


12.64 ±0.08 


7 


16 39 03.30 


-46 42 29.59 


J 


15.99 ±0.07 


14.69 ±0.10 


13.58 ±0.15 


8 


16 39 03.45 


-46 42 27.34 


K, 


16.42 ±0.17 


14.80 ± 0.21 


13.04 ±0.17 


9 


16 39 03.70 


-46 42 27.54 


J 


16.37 ±0.11 


15.00 ±0.11 


13.40 ±0.11 



Astrometric errors are: Ace < O.OP, A5 < 0.1". 



A(v)=10mai 




NTT data 

Giant stars 

Main sequence stars 

Supergiant stars 



1.5 



H - Kshort 



Fig. 6. Near-infrared color-color plot (7 - //) vs (// - K^) de- 
rived using the magnitudes measured by us with the NTT, listed 
in Table 3. The dotted curves show the loci of dwarf, giant 
and supergiant stars, based on the intrinsic colors by Tokunaga 
(2000). The straight dotted lines represent the reddening vectors 
for these stars (see text). 



tainty grows for less intense profiles and largely increases for 
the CO lines, for which no profile fitting is possible as we are 
dealing with unresolved absorption bands. 

Besides identifying hydrogen emission and CO absorption 
features that are typical of pre-main-sequence (PMS) interac- 
tion, we were also able to identify a few absorption lines, which 
can be attributed to the PMS star itself. As shown in Fig. 7, Mg I 
(1.71 jum) and Na I (2.20 /im) seem to be present in our NIR 
spectrum, although somewhat noisy, with equivalent widths of 
~ 3A and > 2A, respectively. Notwithstanding the presence of 
spectral veiling from the circumstellar disk, we can perform a 
first-order approximation of the stellar spectral type, making 
use of the atlas of Rayner et al. (2009) and its discussion in Bik 
et al. (2010). From the cool star atlas, we find that the present 
combination of equivalent widths can only be found in a late 
(K5-M0) dwarf star. 

Source #8. The position of the MSX source G337.9947+0.0770 
(Sect. 2) corresponds to NTT source #8. Mottram and co- 
workers (2010) have measured the 70-;um flux toward this MSX 
source, from the Spitzer-MIPSGal data, also noting that IRAS 



Table 4. Near-IR spectral lines detected in our NTT spectrum 
of source #5. Measured and laboratory wavelengths are listed 
in angstroms in cols. 2 & 3, and their measured radial velocities 
(vlsr), are shown in col. 4. Radial velocities derived for CO lines 
are indicated in brackets as our spectral resolution does not al- 
low us to identify individual lines that constitute the absorption 
band. We have therefore manually measured the wavelength of 
the edge of the absorption profile and used it as a reference value 
for the band head. 



Line ID 




A 


^0 


l^'LSR 






(A) 


(A) 


(kms-') 


Pa (5 




10045.0 


10052 


-220 


Pa/3 




12814.5 


12822 


-186 


'2C0(3- 


-0) 


15578 


15582 


(-88) 


'2CO(4- 


-1) 


15776 


15780 


(-87) 


'^CO (5 - 


-2) 


15978.2 


15982 


(-82) 


'^CO (8 - 


-5) 


16609.2 


16610 


(-29) 


Br 11-4 




16800.5: 


16811 


-198 


Mgl 




17113.3 


171200: 


_ 


Br 10-4 




17369.3 


17367 


-127 


Br 7 




21651.4 


21661 


-137 


NaT 




22050: 


22062 


- 


Cal 




22620: 


22614 


- 


'^CO (2 - 


-0) 


22928.6 


22935 


(-95) 


'^CO (3 - 


-1) 


23222.4 


23227 


(-71) 


"C0(2- 


-0) 


23448 


23439.1 


(-125) 


'^CO (4 - 


-2) 


23516.8 


23535 


(-243) 


"CO (3 - 


-1) 


23728.2 


23739 


(-148) 


'^CO (5 - 


-3) 


23819 


23829 


(-156) 


'^CO (6 ^ 


-4) 


24130.3 


24142 


(-156) 



60 and 100-//m fluxes can be considered upper limit fluxes of the 
same source. We use these mid-infrared fluxes, together with 
our near-infrared photometry, as inputs for the SED builder 
and fitting tool by Robitaille et al. (2007fl The authors have 
built a grid of precomputed radiative transfer models that use a 
fast;\f^ minimization algorithm. If A^ is the number of data points 
(excluding the upper limits), we select the two models for which 
X^ - X^ ■ < 2A^. In Figure 8 we plot the best model that ful- 
fills the last condition. The SED has been built with the NTT 
JHK^ fluxes, the MSX fluxes, the 70//m-Spitzer-MIPS flux, and 
the IRAS 60 and 100//m fluxes as upper limits. We have set a 
distance range from 6 to 10 kpc (see Sect. 4), and an interstellar 
visual absorption range from 5 to 12 (Sect. 3.3). 



Available on-line at |http://caravan.astro.wisc.edu/protostars/| 
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1-55 16 165 

wavelength [^m] 
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Fig. 7. Near-infrared spectrum of the central region of 
IRAS 16353-4636, acquired with the NTT on 2004 July 10. The 
identified emission lines are indicated, most notably Pa^, Pay, 
'^CO lines and bands, '^CO bands, Bry, and possibly Br 10, and 
Br 11. Upper panel: NIR /-band spectrum. Central panel: NIR 
//-band spectrum. Lower panel: NIR K^-hand spectrum. 



4. Discussion 

The H i-line self-absorption technique is a useful tool when de- 
riving distances (see, for instance, Busfield et al. 2006). In this 
way, we measure the most negative absorption feature of the 
VLA-data spectrum on IRAS 16353-4636 (Fig. 4) at v„eg ~ 
-120 km s '. The Galactic rotation model of Brand & Blitz 
(1993) provides two kinematical distances: 6.8 and 9 kpc. The 




Fig. 8. SED fit, using the web tool as explained in Robitaille et 
al. (2007). The dots represent flux values from NTT source #8 
(this work), MSX6C, 70/im Spitzer-MIPSGal, and 60+lOOyum 
IRAS data. 



presence of absorption at -120 km s"' indicates that the source is 
located at or beyond the sub-central point at a distance of ~8 kpc. 
At such a distance, 1" corresponds to 0.04 pc, and the ~10" size 
of the radio source, to 0.4 pc. These are typical values for the size 
of a molecular cloud in the process of contraction to form stars 
(e.g. Garay et al. 2003). The lack of H i absorption at positive 
velocities implies that the source is not beyond the solar circle, 
which sets an upper limit of <15 kpc to its distance. 

The '^CO line results from Urquhart et al. (2007) confirm 
that there is molecular gas with v - Vneg, at least related to NTT 
source #8. The Mopra profile shows five emission lines, with 
LSR-velocities from -122.7 to -38.6 km s"'. Table 5 lists the 
components, widths, and kinematic distances derived from them 
(Urquhart et al. 2007). The central velocity of the strongest line 
is in very good agreement with the extreme H i absorption fea- 
ture at -120 km s '. The rest of the detected '^CO lines can 
be explained, for example, if additional gas clouds are present 
in the direction of IRAS 16353-4636. Observations with better 
angular resolution can help to clarify this issue. 

The velocities measured from the near-IR emission and ab- 
sorption lines in Table 4 are generally consistent with the large 
negative velocities measured in '^CO. Note in particular that the 
strongest emission line, Bry, peaks at an lsr velocity of -137 km 
s ' (i.e. within 15 km s ' the most negative '^CO component). 



Table 5. Parameters derived from the '^CO profile pointing at 
source #8 (20" x 20" angular resolution, Urquhart et al. 2007). 



vlsr 


FWHM 


Tdv 


near-d 


far-rf 


(km/s) 


(knVs) 


(K km/s) 


(kpc) 


(kpc) 


-122.7 


3.3 


29.6 


6.9+0.74 


8.9±0.76 


-70.3 


6.3 


10.3 


4.6±0.44 


11.1 ±0.44 


-57.4 


6.2 


9.8 


4.1 ±0.48 


11.7±0.48 


-47.3 


6.4 


20.1 


3.5±0.54 


12.2±0.54 


-38.6 


7.8 


20.6 


3.0+0.06 


12.8±0.60 



The bolometric luminosity of IRAS 16353-4636, derived 
using the IRAS flux values, is 



--IRAS 



D 



8 kpc 



6.1xlO^Lo, 



Benaglia et al.: The star forming region IRAS 16353-4636 



which is a typical value for massive YSOs (Garay et al. 2003). 

On the other hand, if we assume that the continuum emission 
from NTT source #8 has an optically-thin free-free component 
with flux in the cm wavelengths of the order of 10 mJy (the to- 
tal flux of the entire radio source at 6cm is ~ 80 mJy), and an 
electron temperature of 10"* K for the ionized gas, we find that 
an ionizing photon flux of A^i ~ 5 x 10**^ s"' (Carpenter et al. 
1990) is required to maintain the ionization. This ionizing pho- 
ton flux can be provided by an 09.5 ZAMS star (Panagia 1973). 
The luminosity of such a star is 3.8 x 10"* L©, within a factor or 
two of the bolometric luminosity estimated from the IRAS fluxes. 
Although the FIR colors could be compatible with an ultracom- 
pact H II region, the source does not appear in the CS survey by 
Bronfman et al. (1996). 

The good fit of the SED model (Fig. 8) of the infrared spec- 
trum of source #8, together with its location in the NIR C-C 
diagram (Fig. 6), and the agreement with the coordinates of the 
IRAS and MSX sources in the region, make it the most promis- 
ing candidate for a MYSO in this cluster. The best fit for an 
object at ~ 8 kpc, with Ay ~ 10, is achieved for a model source 
of ~ 23 solar masses, ~90 solar radius, ~ 9000 K effective tem- 
perature, and 5000 yr of age. 

5. Conclusions 

From the multi-frequency study presented here on IRAS 
16353-4636, we have found strong evidence that supports its 
nature as a star-forming region. The H i absorption spectrum in- 
dicates that the source is located at or beyond the sub-central 
point at a distance of ~8 kpc. At such a distance, the size of the 
radio source agrees with typical values for the size of a molecu- 
lar cloud in the process of contraction to form stars. The region 
harbors, at least, two pre-main sequence stars or YSOs (NTT 
sources #5 and #8), plus three YSO candidates (NTT sources #6, 
#7, and #9). Our photometry shows that these five point sources 
have infrared excess. NIR spectra obtained from one of these 
(#5) indicate that it is a low-mass pre-main sequence star, whilst 
the combination of near- and mid-infrared data suggests that 
another one (#8) is a massive young stellar object. 

The broad H i absorption feature and multiple CO emission 
lines seen at various velocities along the line of sight towards 
this region suggest the presence of multiple cloud components. 

The radio results reveal that the region is complex, with fine- 
scale structure seen at the highest frequencies. The results of a 
spectral index analysis at mm wavelengths are indicative of dif- 
ferent physical conditions varying according to position. From 
1.4 to 19.6 GHz, the radio emission presents a peak, character- 
ized with a negative spectral index. This can be explained if the 
peak represents the terminal point of an outflow. 

Obtaining NIR spectra for the remaining sources in this re- 
gion will allow us to improve our knowledge of the whole star 
forming region, as high-quality spectra provide both physical in- 
sight into their sources and kinematical information necessary 
to disentangle projection effects. Higher angular resolution and 
better sensitivity observations at low radio frequencies to search 
for polarized emission, combined with molecular line and con- 
tinuum studies, are needed to deepen our investigation of this 
interesting source. 
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